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ABSTRACT 

The well studied blazar Markarian 421 (Mrk 421, 2;=0.031) was the subjeet of an 
intensive multi-wavelength campaign when it flared in 2013 April. The recorded X-ray 
and very high energy (VHE, E>100 GeV) y-ray fluxes are the highest ever measured 
from this object. At the peak of the activity, it was monitored by the hard X-ray focus¬ 
ing telescope Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift X-Ray Tele¬ 
scope (XRT). In this work, we present a detailed variability analysis of NuSTAR and 
Swift-XKT observations of Mrk 421 during this flaring episode. We obtained the short¬ 
est flux doubling time of 14.01±5.03 minutes, which is the shortest hard X-ray (3—79 
keV) variability ever recorded from Mrk 421 and is on the order of the light crossing 
time of the black hole’s event horizon. A pattern of extremely fast variability events 
superposed on slowly varying flares is found in most of the NuSTAR observations. We 
suggest that these peculiar variability patterns may be explained by magnetic energy 
dissipation and reconnection in a fast moving compact emission region within the jet. 
Based on the fast variability, we derive a lower limit on the magnetic field strength of 
B > 0.735]”^^^ G, where is the Doppler factor in units of 10, and z/ig is the 
characteristic X-ray synchrotron frequency in units of 10^® Hz. 

Subject headings: galaxies: active — X-rays: galaxies — quasars: individual (Mrk421) — 
galaxies: jets 
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1. Introduction 

Blazars are a special class of radio-loud active galactic nuclei (AGN) whose observed 
broadband spectra are dominated by highly variable, non thermal, and Doppler boosted radiation 
from powerful relativistic jets. An exhaustive and detailed search for blazar variability at different 
wavelengths is necessary to understand not only the size and/or location of the emission region, 
but also the involved particle acceleration mechanisms and radiative processes. In this aspect, 
observations of extremely fast variability at very high energies (VHE, E> 100 GeV) from some BE 
Eac objects (e.g., Aharonian et al. 2007) have raised questions on the radiative models proposed 
to explain blazar emissions. However, the lack of sensitive hard X-ray instruments prevented us 
to understand the connection between the observed variability at VHE, corresponding to the tail 
of the synchrotron self Compton (SSC) spectrum, with that at hard X-rays (synchrotron tail), 
under the leptonic emission scenario. Thanks to the hard X-ray mission Nuclear Spectroscopic 
Telescope Array (NuSTAR, Harrison et al. 2013), such valuable information at hard X-ray energies 
are now available and using them one can get deeper insights about the physics of blazar radiation 
zones. 

Markarian 421 (Mrk 421, z = 0.031) is a BE Eac object which has been studied extensively 
over a broad spectral range starting from radio to VHE 7 -rays (see e.g., Eossati et al. 2008; 
Abdo et al. 2011; Shukla et al. 2012; Aleksic et al. 2015a,b). The lack of emission lines and a 
thermal component in its broadband spectral energy distribution (SED) suggests that the emission 
from the jet to be dominant rather than other external sources like the accretion disk, and the 
broad line region. Consistently, the high energy emission is often explained successfully by an 
SSC process without invoking any additional radiative component (e.g., Aleksic et al. 2015b). In 
addition, the extension of the synchrotron spectrum up to X-rays reflects an efficient acceleration 
mechanism that produces a particle spectrum extending up to extremely relativistic energies. 
Accordingly, the SSC spectrum also extends to VHE 7 -rays and hence, Mrk 421 is known to be a 
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strong TeV emitter (Punch et al. 1992). It exhibits a flat radio spectrum, optical polarization, and 
large amplitude variability throughout the electromagnetic spectrum (e.g., Aleksic et al. 2015a; 
Hovatta et al. 2015). In particular, extremely fast VHE outbursts were detected where the doubling 
time of the flare events were found to be <15 minutes (Gaidos et al. 1996). A positive correlation 
between X-rays and VHE radiation is also reported (Aleksic et al. 2015a), thus suggesting these 
emissions arise from the same region. Using the Whipple observatory {E > 400 GeV) and 
All-Sky Monitor (ASM, 2—10 keV) onboard Rossi X-ray Timing Explorer (RXTE) data, a long 
term study of Mrk 421 has been performed by Acciari et al. (2014) and they also report a positive 
correlation between VHE and X-ray fluxes. 

In 2013 April, Mrk 421 underwent a giant X-ray flare and was extensively monitored by 
both space and ground based observational facilities (e.g., Cortina & Holder 2013; Paneque et al. 
2013), including NuSTAR. In this paper, using publicly available NuSTAR and Swift X-Ray 
Telescope (XRT) data, we perform a detailed study of the X-ray emission (0.3—79 keV) from 
Mrk 421 covering the period of this exceptional outburst. A major emphasis is given on searching 
for the fastest variations seen in this energy regime and also on understanding the patterns of hard 
X-ray variability. We use a ACDM cosmology with the Hubble constant Hq = 71 km s“^ Mpc~^, 
Vtm = 0.27, and Ha = 0.73. 


2. Observations 

2.1. NuSTAR 

NuSTAR is a hard X-ray focusing satellite consisting of two co-aligned focal plane modules 
referred as EPMA and EPMB* (Harrison et al. 2013). Mrk 421 was observed by NuSTAR first 


'https://heasarc.gsfc.nasa.gov/docs/nustar/ 



-5- 


in 2012 for pointing calibrations and later in 2013 as part of a coordinated multi-wavelength 
campaign (Balokovic et al. 2013). During the peak of the 2013 April outburst, NuSTAR monitored 
Mrk 421 many times between April 10 to April 20 (MJD 56392—56402) and we define this as 
the period of interest. The NuSTAR data are analyzed with the package NuSTARDAS (v. 1.4.1). 
Calibration and cleaning of the event files are done with the task nupipeline and using 
NuSTAR CALDB updated on 2015 January 23. In the energy range of 3—79 keV, the source 
light curves and spectra are extracted using nuproduct s, from a circular region of 30" radius 
centered at the position of Mrk 421, whereas the background region is selected as a circle of 70" 
radius from a nearby region free from source contamination. To generate light curves, both FPMA 
and FPMB count rates are summed and background subtracted. To develop the light curves 
corresponding to one data point per observation ID, we select the bin size as the total duration of 
the observation run, whereas finer binned light curves are extracted using time bins as short as 5 
minutes. In principle, one can adopt even shorter time binning, but choosing extremely short bins 
may result in larger flux uncertainties and a poorly constrained flux doubling time. On the other 
hand, a longer bins will wash out short timescale features. With this in mind, we find that 5 minute 
binning is optimal. Source spectra are binned to have at least 20 counts per bin and spectral fitting 
is performed in XSPEC (Amaud 1996). 


2.2. Swift-XRT 

Swi/t-XRT (Burrows et al. 2005) data, covering the period of the outburst, have been analyzed 
using the publicly available “Sw///-XRT data product generator” available at the University of 
Leicester website^. The details of the methodology adopted can be found in Evans et al. (2007, 
2009). We extract the XRT light curves using the bin size equal to the total exposure of a particular 

^http ://w w w. s wift. ac. uk/u ser_obj ects/ 



- 6 - 


observation, and also with 5 minute time binning. Moreover, light eurves are also generated in 
three different energy bands with 5 minutes binning: 0.3—10 keV, 0.3—1.5 keV, and 1.5—10 keV. 
In both NuSTAR and Swift-XKT analysis, we rejeet all the bins with Fx/AFx < 3, where AFx is 
the assoeiated error in the X-ray flux Fx- 


3. Results 

The NuSTAR and Swi/t-XRT light eurves of Mrk 421, eovering the period of high aetivity, 
are presented in Figure 1. In this plot, each data point corresponds to one observation ID. As 
can be seen, two prominent flares are visible, one at around MID 56395 and another at around 
MID 56397. Further, the good photon statistics during the outburst permits us to generate light 
curves using shorter binning as small as 5 minutes. Generation of such shorter time binned light 
curves is also useful in searching for faster variability and the possible existence of patterns in the 
flux variations. The resultant plots are shown in the top panel of Figure 2 and 3 for NuSTAR and 
Swi/t-XRT data, respectively. Moreover, we also generate the light curves in two different energy 
bands, both for NuSTAR (3-10 keV and 10-79 keV) and Swift-XKY (0.3-1.5 keV and 1.5-10 
keV), and they are shown in the middle panels of Figure 2 and 3. These observations indicate the 
presence of intra-day variability. Moreover, as can be seen, multiple episodes of flaring activities 
are observed both in soft and hard X-ray bands. For better visualization of the patterns of flux 
variations, 5 minute binned NuSTAR light curves are also presented in Figure 4. In this plot, 
each panel represents one individual NuSTAR pointing. As can be seen in Figure 4, and also in 
Figure 2, two distinct patterns are visible, a slowly varying flare and on the top of that extremely 
fast flux variations. Though there are also hints for similar behavior in the XRT light curves, a 
strong claim cannot be made due to the short exposure of the XRT observations. 


We calculate the normalized rms amplitude of variability (F^ar, Vaughan et al. 2003) to 
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characterize the flux variations. It is defined as follows: 


/ 


(F)2 


( 1 ) 


where (F) is the mean flux, S'^ is the the sample variance, and is the mean-square value of 
the uncertainties. The error in F^ar is computed as follows (Poutanen et al. 2008; Aleksic et al. 
2015a), 



( 2 ) 


Fvar is calculated for all the light curves shown in Figure 2 and 3. For the NuSTAR observations 
overall in the 3—79 keV energy range, the Fvar is found to be 0.790 ± 0.001, whereas, for the 
XRT data (0.3—10 keV) it is 0.599 ± 0.007. Due to the long exposure of each NuSTAR pointing, 
we are able to derive F^ar for individual NuSTAR observation id also. The results are presented in 
Table 1. It is important to note here that the F^ar for 10—79 keV is found to be higher than that 
for 3—10 keV in almost all the observations. Similar behavior is noticed for the Sw///-XRT data 
where a Fvar of 0.565 ± 0.006 and 0.693 ± 0.006 is obtained for the 0.3—1.5 keV and 1.5—10 
keV energy ranges, respectively. Further, the fine binned light curves presented in the top panel 
of Figure 2 and 3 are scanned to find the shortest flux doubling/halving time using the following 
formula 


F{t) = 


(3) 


where r is the characteristic doubling/halving timescale and F(to) and F{t) are the fluxes at time 
to and t, respectively. The conditions that the difference in fluxes at the epochs t and to is at least 
3a significant is also met (see e.g., Foschini et al. 2011). The shortest flux doubling time for XRT 
data is 64.14 ± 13.78 minutes with ~5cr significance. The absence of minute scale variability 
(< 15 minutes) in fine binned XRT light curves of Mrk 421 was earlier reported by Pryal et al. 
(2015). However, the fastest flux doubling time (fvar) for NuSTAR observations is found to be 
14.01 ±5.03 minutes. If a more conservative 5cr significance is considered, then the shortest flux 
doubling time is 28.44±3.76 minutes. This is the shortest hard X-ray variability ever detected 
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from Mrk 421 and is on the order of the light crossing time of the black hole’s event horizon (see 
Section 4). The parameters associated with this analysis are given in Table 1. 


In the bottom panels of Figure 2 and 3, we show the temporal variations of the hardness ratio 
(HR). It is calculated using the following equation 


HR 


-^hard .^soft 

-^hard T -^soft 


(4) 


where Fsoft and Fhard are soft (0.3—1.5 keV for XRT and 3—10 keV for NuSTAR) and hard 
(1.5—10 keV for XRT and 10—79 keV for NuSTAR) X-ray fluxes, respectively. A ‘harder when 
brighter’ behavior is evident for both NuSTAR and XRT light curves, especially at the peak of the 
flare around MJD 56397. 


To determine the spectral behavior of the source, we fit the NuSTAR spectra with a log 
parabola model (see e.g., Massaro et al. 2004), while keeping the pivot energy fixed to 10 keV. 
The results of the spectral fitting are presented in Table 2. Moreover, the variations of both 
spectral parameters, i.e. the photon index («) at the pivot energy and the curvature index (/3), as 
a function of brightness are also shown in Figure 4. The spectra are found to be more curved 
during higher flux states, whereas there is a clear trend of ‘hardening when brightening’, thus 
supporting the behavior seen in the hardness ratio plots. We note that the joint Sw(//-XRT and 
NuSTAR spectral fitting, for the same period, has recently been performed by Sinha et al. (2015) 
and thus it is not presented here. However, the spectral behavior of Mrk 421 observed from joint 
XKT-NuSTAR spectral analysis, as done by Sinha et al. (2015), is similar to that obtained by us. 


4. Discussion and Conclusions 

The blazar Mrk 421 is known to exhibit fast variability at all wavelengths, especially at 
X-rays and VHE 7 -rays (Gaidos et al. 1996; Cui 2004). Though a positive correlation between 
these two energy bands is frequently found (e.g. Aleksic et al. 2015a), simultaneous hard X-ray 
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and VHE observations were lacking during earlier measurements. The simultaneity becomes more 
important considering that, in leptonic models, both hard X-ray and VHE photons are expected 
to be produced by the same population of high energy electrons. The recent X-ray outburst of 
Mrk 421 was contemporaneously monitored by NuSTAR and ground based Cherenkov telescopes, 
thus providing an excellent opportunity to constrain the radiative processes in a way that was not 
possible before. 

Using RXTE observations of Mrk 421, Cui (2004) reported the presence of minute scale 
X-ray variability, however, they could not quantify the parameters due to gaps in the data. 
Recently, Pryal et al. (2015) have searched for fast X-ray variability (<15 minutes) among a 
sample of AGN monitored by the 5w///-XRT, but were unsuccessful. In the energy range of 2—10 
keV, the RXTE light curves of Mrk 421, covering the entire duration of RXTE monitoring, are 
publicly available^ and the details of the data reduction procedure are provided in Rivers et al. 
(2013). Using these results, we calculate the shortest flux doubling/halving time and found it to 
be 1.38 ± 0.37 hr. Therefore, to our knowledge, this is the first time that a statistically significant 
hard X-ray flux variability, as small as ~14 min, has been detected from Mrk 421. 

The shortest hard X-ray variability time estimated in this work is ~14 min, which is similar 
to that observed in the VHE band by Gaidos et al. (1996). Interestingly, during the 2013 April 
outburst, Mrk 421 seems to show fast variability also at VHE 7 -rays (Cortina & Holder 2013). 
This observation, therefore, suggests a co-spatial origin of the X-ray and 7 -ray flares. If the black 
hole mass of Mrk 421 is taken as 1.9 x 10® Mq (Barth et al. 2003), the observed hard X-ray 
variability timescale is approximately identical to the light-crossing timescale across the black 
hole’s event horizon (feH ~ ^g/c = GM/ c® ~ 15 min) which is the shortest expected variability 
timescale of emission powered by accretion onto the black hole. The variability timescales 
estimated from NuSTAR observations are, thus, difficult to explain using the conventional blazar 


®http://cass.ucsd.edu/~rxteagn/ 
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radiation models. 

The detection of extremely fast variability from several blazars seriously challenges 
the commonly accepted single-zone jet models for blazar emission (e.g., Albert et al. 2007; 
Aleksic et al. 201 1). In the framework of such models, the Doppler factor of the compact emitting 
region has to be very high (>50) in order to avoid the severe pair production of TeV photons with 
the synchrotron radiation and in some cases (although not in the case of the X-ray variability of 
Mrk 421 presented here) also to satisfy the condition fvar ^ ^bh (e.g., Begelman et al. 2008). 
However, interferometric observations of superluminal radio knots suggest lower values of the 
Doppler factor (Lister et al. 2009). This apparent contradiction can be avoided by arguing that 
radio and hard X-ray emissions come from different emission regions and the jet is decelerated at 
sub-pc scales (e.g., Levinson 2007) after the production of the TeV emission. On the other hand, 
several alternative models have been proposed to explain such fast variability (of VHE radiation, 
in particular). Ghisellini & Tavecchio (2008) have invoked the localized magneto-centrifugal 
acceleration of beams of electrons to explain the fast TeV variability of PKS 2155—304 and 
Mrk 501, however, this ‘needle’ model predicts little or no variability in X-rays. The high 
activities seen in both VHE and X-rays from Mrk 421 during 2013 April outburst disfavors 
this hypothesis. Giannios et al. (2009) proposed a ‘jet-in-a-jet’ model, in which the concept 
of magnetic reconnection is used to explain the observed fast variability. This model not only 
reproduces the extremely fast TeV variations, but also predicts the observations of fast X-ray 
flares. The observed extremely fast hard X-ray flux variations along with the hint for high 
flux activity at TeV energies (Cortina & Holder 2013) strengthens the hypothesis that magnetic 
reconnection is a possible origin of the 2013 April flare of Mrk 421. Moreover, the model of 
Giannios et al. (2009) also predicts the presence of a slowly varying flare due to tearing of a 
large reconnection region. This leads to the ejection of several individual relativistic plasmoids 
which are thought to be responsible for fast variations. As can be seen in Eigure 2 (and also 
in Eigure 4), we do see slowly varying patterns underlying more rapid, short-term flares. This 
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provides further support for the hypothesis of magnetic dissipation. It is interesting to note here 
that the above mentioned models (see also Narayan & Piran 2012) have a common assumption 
of a small emission region moving much faster than the surrounding jet medium. A rapid flare 
can be observed, thus, by shorter light crossing timescale along with stronger beaming effects. 
On a completely different note, Zacharias (2014, and references therein) has invoked the time 
dependent particle injection with non-linear SSC cooling to explain the fast variability seen in the 
blazar light curves. 

Though the fundamental causes of the origin of the 2013 April outburst of Mrk 421 are 
uncertain, a few model-independent parameter estimates can be derived merely based on the 
assumption of a synchrotron origin of the hard X-ray emission from Mrk 421. If the radiation 
output of the dominant electron population is primarily by synchrotron emission, electrons of 
energy lose energy on an observed timescale of tc = ([1 -f z\/5) {Qtt nieC^) / {caT 

where <5 = 10 (5i is the bulk Doppler factor, B = 1 Bq Gauss is the magnetic field, and ctt is the 
Thomson cross section. The electron Lorentz factor can be associated with a characteristic X-ray 
frequency in the NuSTAR energy range, z/gy = 10^® z/ig Hz = 4.2 x 10® {5/[l + z\) Bq 7 ^ Hz. 
Combining these two identities and requiring that the synchrotron cooling timescale of electrons 
radiating in the NuSTAR regime has to be shorter than or equal to the observed minimum 
variability timescale, we find 


H > 0.73 5”^/^ z/if G 


(5) 


1 /3 

Thus, even for 5 ~ 30, the inferred magnetic field of i? > 0.35 G is higher than the values of 


H < 0.1 G typically inferred from SED modeling of high frequency peaked BE Eacs such as Mrk 
421. Assuming values oi B = 0.4 G and 5 = 30, consistent with the above estimates, electrons 
radiating near the high-energy end of the NuSTAR range, have Eorentz factors of 



( 6 ) 


and can produce 7 -rays by Compton scattering in the Thomson regime up to photon energies of 
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,max 


((5/[l + z]) 'yrUeC^ or 



1/2 


B 


1/2 


ET,max ~ 6 


,max 


0.4 G 



(7) 


by scattering target photons of energy Et ~ 38 (5/30)^/^ (i?/0.4G)^/^ eV, i.e., UV - 
soft X-ray photons. Henee, the same population of ultrarelativistic eleetrons ean plausibly be 
responsible for both hard X-ray synehrotron and Compton VHE 7 -ray emission, varying on 
eomparable timeseales, thus providing strong support for a leptonie (plausibly synehrotron self 
Compton) co-spatial origin of the X-ray and VHE 7 -ray emission. 

The extremely fast variability seen at hard X-rays suggests the impulsive injeetion 
(aeeeleration) of eleetrons of the highest energies as the likely eause of the flux variations, sinee 
the highest energy eleetrons have the shortest eooling timeseales. The injeetion of highly energetie 
partieles is expeeted to eause not only a flux enhaneement but also a speetral hardening, whieh is 
seen. However, what eauses the injeetion of the highest energy eleetrons and/or what ean energize 
the particles remains unclear. 
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Table 1: Variability Characteristics of Mrk 421 for 5 minutes binned NuSTAR Observations shown 
in Figure 2. Col.[l] and [2]: NuSTAR observation date and observation ID; Col.[3], [4], [5]: nor¬ 
malized rms variability amplitude for 3—79 keV, 3—10 keV, and 10—79 keV respectively; Col.[6] 
and [7]: the shortest flux doubling/halving time in minutes and its significance obtained for 3—79 
keV data; Col.[8]: R: rising flare D: decaying flare. 


Date 

[1] 

OBSID 

[2] 

Tp3 — 79 keV 
var 

[3] 

p3 — 10 keV 

[4] 

plO —79 keV 
-^var 

[5] 

kl 

[6] 

Signif. 

[7] 

R/D 

[8] 

2013-04-10 

60002023024 

0.150 ±0.005 

0.149 ±0.007 

0.146 ±0.024 

159.56 ±40.13 

4.17 

D 

2013-04-11 

60002023025 

0.604 ± 0.002 

0.595 ± 0.002 

0.642 ± 0.006 

14.01 ±5.03 

3.15 

R 

2013-04-12 

60002023027 

0.140 ±0.002 

0.137 ±0.003 

0.185 ±0.008 

47.97 ±4.13 

11.33 

R 

2013-04-13 

60002023029 

0.254 ± 0.002 

0.223 ± 0.003 

0.219 ±0.007 

- 

- 

- 

2013-04-14 

60002023031 

0.320 ± 0.001 

0.317 ±0.002 

0.366 ± 0.003 

37.87 ± 4.96 

7.63 

R 

2013-04-15 

60002023033 

0.184 ±0.008 

0.186 ±0.002 

0.221 ± 0.006 

36.46 ± 9.77 

3.73 

D 

2013-04-16 

60002023035 

0.400 ± 0.002 

0.398 ± 0.002 

0.451 ±0.005 

28.44 ± 3.76 

7.57 

R 

2013-04-18 

60002023037 

0.186 ±0.005 

0.182 ±0.005 

0.227 ± 0.015 

20.26 ±5.71 

3.55 

R 

2013-04-19 

60002023039 

0.123 ±0.006 

0.121 ±0.006 

0.157 ±0.018 

- 

- 

- 
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Table 2: Summary of the NuSTAR Data Analysis 


OBS ID^ 

Obs. date^ 

Exp.^ 


/3" 

7^3-79 keV® 

Normalization^ 

Stat.8 

60002023024 

56392.89 

5.76 

3.011 ±0.022 

0.326 ± 0.057 

5.883 ± 0.066 

1.250 ±0.013 

634.40/624 

60002023025 

56393.04 

57.51 

2.725 ± 0.005 

0.298 ±0.012 

10.460 ± 0.028 

2.383 ± 0.006 

1799.28/1412 

60002023027 

56394.86 

7.63 

2.735 ± 0.009 

0.388 ± 0.024 

22.380 ±0.117 

5.217 ±0.024 

1084.49/1014 

60002023029 

56395.90 

16.51 

2.908 ±0.011 

0.338 ± 0.029 

7.790 ± 0.055 

1.716 ±0.009 

1003.08/911 

60002023031 

56396.90 

15.61 

2.390 ± 0.005 

0.360 ±0.013 

28.605 ±0.106 

6.819 ±0.020 

1715.36/1424 

60002023033 

56397.92 

17.28 

2.672 ± 0.009 

0.283 ± 0.024 

9.088 ± 0.047 

2.080 ±0.010 

1036.17/1012 

60002023035 

56398.93 

20.28 

2.466 ± 0.007 

0.287 ±0.019 

11.021 ±0.056 

2.570 ±0.011 

1231.30/1188 

60002023037 

56400.01 

17.80 

2.966 ± 0.027 

0.290 ± 0.068 

1.384 ±0.015 

0.296 ± 0.004 

526.60/557 

60002023039 

56401.02 

15.96 

3.031 ±0.031 

0.179 ±0.076 

1.279 ±0.014 

0.259 ± 0.004 

515.62/523 


^NuSTAR observation id. 

^Date of observation, in MJD. 

^Net exposure, in ksec. 

^Photon index at pivot energy, in the log parabola model. 

^Curvature index, in the log parabola model. 

^Energy flux in 3—79 keV energy band and in units of 10“^° erg cm“^ s“^. 
^Normalization in units of 10“^ ph cm“^ keV“^. 

^Statistical parameters: x^ldof. 
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Fig. 1.— Flux history of Mrk 421 during the period of high aetivity, as observed from NuS¬ 
TAR (3—79 keV, top panel) and Sw(//-XRT (0.3—10 keV, bottom panel). In both the panels, eaeh 
data point corresponds to one observation ID. 
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Fig. 2.— NuSTAR light curves of Mrk 421 in the energy range of 3—79 keV (top), 3—10 keV 
(seeond from top), and 10—79 keV (second from bottom). Bottom panel refers to the variation of 
hardness ratio (dehned as, HR = ). The fluxes are in units of eounts s“^ and the adopted 

time binning is 5 minutes. 
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Fig. 3.— Swz//-XRT observations of Mrk 421 during the period of outburst. The seleeted time 
binning is 5 minutes and the fluxes are in units of eounts s“^. The variation of the hardness ratio 
(HR = ) is shown in the bottom panel. 
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Fig. 4.— NuSTAR light curves of Mrk 421, same as in the top panel of Figure 2, but zoomed 
to show the pattern of variations . Eaeh panel represents the NuSTAR observation taken on that 
specihe day. Other information are same as in Figure 2. 


















Photon index (a) 


-23- 



Fig. 5.— Variations of log parabolic spectral parameters as a funetion of 3—79 keV energy flux. 
Left y-axis represents the photon index at the pivot energy whereas right y-axis eorresponds to 
eurvature index. 


Curvature index [( 5 ) 






















